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Abstract 

The  constitutive  equations  for  materials  with 
memory  developed  by  Green  and  Rivlin  are  discussed, 
with  particular  emphasis  on  the  physical  assumptions 
regarding  material  behavior  which  are  implied  by  the 
mathematical  assumptions  in  the  theory. 
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1.  Introduction 

The  ciiarsctcristr r  property  of  viscoelastic  soi-«.s  which 
distinguishes  them  from  perfectly  elastic  eiids  is  the  fact 
that,  it  they  are  subjocteu  to  a  deform  tion  which  varic:  with 
time,  the  stress  bk  ^urtc  at  time  t  ,  t>av.  depends  not  only 
on  the  instantaneous  value  of  th-3  defot nation  gradients,  but 
also  on  the  whole  previous  his„oiy  cf  the  deformation  g~r  cl¬ 
ients.  In  a  series  of  papers  Green  and  Riv» \n  (195/,  196u) 
and  Green,  Spencer  and  Rivlin  (1959)  have  developed  constitu¬ 
tive  equations  for  such  materials,  in  which  the  stress  is  ex¬ 
pressed  in  terms  of  the  deformation  in  the  form  of  series  of 
multiple  integrals.  It  is  the  object  of  this  paper  to  recapitu¬ 
late  this  development,  with  part  cular  emphasis  on  tnc  physical 
assumptions  regarding  the  ..  »..*•?  f.ich  are  implied  by  the 
mariv  jmatical  as  sump  ions  made  i*»  the  theory.  Such  a  develop- 
•  is  perhaps  timely  in  view  of  the  extensive  attempts  in 
recent  years  to  represent  the  behavior  of  actual  materials  in 
the  form  given  by  the  theory. 

In  §2,  we  start  with  tne  assumption  that  the  Cauchy 
stress  is  a  functional  of  the  deformation  gradients  for  times 
up  to  and  including  t  and  express  the  restrictions  imposed 
on  the  form  of  this  dependence  by  the  consideration  that,  if 
the  body  and  the  force  system  applied  to  it,  including 
inertial  forces,  are  simultaneously  subjected  to  a  time- 
dependent  rigid  rotation,  the  stress  field  at  time  t  is 
correspondingly  rotated  by  the  amount  of  this  rotation  at 
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tiae  t  ,  The  fi ro‘.  ional  dependence  on  the  history  of  the 
d‘j£orp*tion  gra^L.07;  5  i  p  to  and  including  time  t  is  thus 
replaced,  wit]-,  any  jsired  approximation,  by  functional 
dependence  on  the  v  ,!uced  Cauchy  strain  (i.e.,  the  Cauchy 
strain  less  fat  7  i  tensor).  The  nature  of  this  dependence 
can  be  further  rest}  cted  if  the  material  has  some  symmetry, 
b\it  the  explicit  dci.  rmination  of  the  form  of  this  restric¬ 
tion  will  not  bo  <ii"  ossed  here.  Rather  we  shall  be  concerned 
with  the  assumptions  regarding  the  functional  dependence  of  the 
stress  on  the  redeem  Cauchy  strain  which  enable  us  to  obtain 
explicit  represent  -  ; : ons . 

Sir*,  much  of  the  apparent  complexity  of  the  theory,  as 
originally  presented,  stems  from  the  tensor  character  of  the 
equations  involved,  we  first  discuss  the  simple  case  in  which 
only  simple  extension  of  a  rod  of  the  material  under  tensile 
force  is  considered.  The  stress  is  then  considered  to  be  a 
f  .notional  of  the  history  of  the  fractional  extension  of  the 
rod  up  to  and  including  time  t  .  The  assumption  is  made 
that,  for  two  extension  histories,  for  which  the  extensions 
at  each  instant  of  time  are  infinitesimally  different,  the 
corresponding  stresses  at  time  t  are  infinitesimally  differ¬ 
ent.  This  assumption  that  the  dependence  of  the  stress  on 
the  extension  history  is  continuous,  enables  us  to  obtain, 
with  any  desired  approximation,  representations  for  the  stress 
at  time  t  as  the  sum  of  a  series  of  multiple  integrals  in 
the  extension  history.  It  is  shown  that  the  kernels  in  these 
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multiple  integrals  must  necessarily  be  such  as  to  yield  the 
property  that  the  memory  of  the  stress  for  extension  history 
is  a  fading  memory  -  the  stress  forgets  the  extension  history 
in  the  distant  enough  past.  This  results  from  the  continuity 
assumption  and  from  the  assumption  that,  whatever  the  exten¬ 
sion  history,  the  stress  does  not  become  infinite  at  infinite 
time. 

It  is  pointed  out  that  the  basic  assumption  that  the 
stress  at  time  t  is  a  functional  of  the  extension  history 
up  to  and  including  time  t  is  strictly  not  broad  enough  to 
accommodate  the  behavior  of  materials  which  h«ve  internal 
friction,  if  deformations  involving  ins  tan taneo  is  changes  of 
the  velocity  gradient  are  allowed.  In  such  cases  the  stress 
in  the  rod  at  time  t  must  be  assumed  to  depend  not  only  cn 
the  history  of  the  extension  up  to  and  including  t  ,  but 
also  on  the  instantaneous  values  of  the  rate  of  extension  at 
time  t  . 

Finally  in  §4,  the  various  types  of  behavior  considered 
in  63,  in  the  case  of  time-dependent  simple  extension  of  a  rod, 
are  generalized  to  the  case  of  arbitrary  deformations  of  the 
material. 
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2.  Basic  theory 

In  this  section,  the  principles  underlying  the  continuum 
mechanical  theory  of  non-linear  viscoelastic  solids  will  be 
briefly  summarized.  We  start  by  delimiting  in  mathematical 
terms  the  class  of  materials  to  which  the  theory  is  to  apply. 
This  class  is  not  the  most  general  with  which  we  shall  be  con¬ 
cerned  in  this  paper.  However,  it  will  serve  to  illustrate  t.»e 
fundamental  principles  involved. 

The  deformation  is  described  by  specifying  the  vector 
position  x(x)  of  a  generic  particle  of  the  body  at  time  x  , 
with  respect  to  a  fixed  origin,  as  a  function  of  its  vector 
position  X  with  respect  to  the  same  origin  at  some  reference 
time  T  ,  say,  which  is  conveniently  taken  to  be  a  time  at 
which  the  body  is  undeformed.  Thus, 

x(x)  =  x(X,x)  .  (2.1) 

The  displacement  vector  u(x),  defined  by 

u(x)  =  x(x)  -  X  ,  (2.2) 

may  also  be  regarded  as  a  function  of  X  and  x  ,  thus 

u(x)  =  u(X,x)  .  (2.3) 

In  terms  of  the  components  x  (x),  X.  and  u (x) 

p  A  p 

o±  x(x),  X  and  u(x)  respectively,  in  a  rectangular 
cartesian  coordinate  system  x  ,  the  relation  (2.1)  may  be 
written 


xp(x)  ■ 


xp(XA,x) 


(2.4) 
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and  the  relation  (2.3)  may  be  written 


Up(T)  *  “p(XA.T)  • 


(2.5) 


It  is  evident  that  if  the  nine  deformation  grad¬ 
ients  3x (t) / 3X-  are  specified  at  a  point,  then  the  amount 

p  A 

by  which  any  linear  element  at  that  point  is  stretched  is 
determined.  For  brevity,  we  use  the  so-called  comma  notation, 
in  which  A  denotes  the  operator  3/3XA  .  The  nine  defor¬ 
mation  gradients  3x  fx)/3X.  are  then  written  x  .  (x)  . 

p  a  P>-ft 

We  shall  be  concerned  with  materials  in  which  the  stress 
Ojj(t)  at  a  generic  particle,  measured  at  time  t,  which,  for 
brevity,  will  b°  denoted  ,  may  depend  on  the*  history  of 

the  deformation  gradients  x  .  (t)  at  that  particle  for  all 

p  ,  A 

times  from  x  =  -  to  t  .  Thus,  unlike  the  situation 
which  exists  in  an  elastic  material,  in  which  the  stress  at 
time  t  is  determined  uniquely  by  the  deformation  gradients 
at  time  t  ,  in  a  viscoelastic  material  the  stress  at  time  t 
"remembers"  the  deformation  gradients  at  all  previous  times 
and  the  material  is  said  to  possess  "memory".  We  can  express 
this  physical  idea  in  a  mathematical  statement  -  the 
constitutive  assumption  for  the  material  -  that  the  stress 
cfjj  at  time  t  is  a  tensor- valued  functional  of  the  deforma¬ 
tion  gradient  history  x_  .(x)  for  x  =  -  to  t,  thus 

P  9* 


(2.6) 


T=  -  oo 


This  means  simply  that  if  the  dependence  of  x_  a(t)  on  x 

pf  A 

is  known  for  a  generic  particle,  then  the  stress  at  that 
particle,  at  time  t,  is  determined. 

We  recall  that,  in  classical  elasticity  theory,  if  we 
make  the  constitutive  assumption  that  the  stress  at  time  t 
depends  on  the  values  of  the  nine  displacement  gradients 
3u^/3xj  at  time  t,  it  is  shown  that  it  must  in  fact  do  so 
through  the  six  infinitesimal  strain  comporents  e^^ 
defined  by 


The  question  arises  -  are  there  any  corresponding  restric¬ 
tions  which,  in  the  case  of  the  constitutive  assumption 
(2.6)  for  a  non-linear  viscoelastic  material,  can  be  placed 
on  the  manner  in  which  o^j  depends  on  the  deformation  grad¬ 
ient  history?  That  such  restrictions  do,  in  fact,  exist 
follows  from  a  consideration  somewhat  similar  to  Oat  used  in 
classical  elasticity  theory.  We  consider  that  the  body,  to¬ 
gether  with  the  force  system  associated  with  it, including  iner¬ 
tial  forces,  is  subjected  to  a  time -dependent  rigid  rotation.  The 
stress  field  at  time  t,  referred  to  the  system  x,  differs  from 
that  which  obtains  ip,  the  absence  of  the  superposed  rotation, only 
to  the  extent  that  it  undergoes  a  rotation,  the  amount  of  which 
is  that  of  the  superposed  rotation  at  time  t.  From  this 
consideration,  it  follows,  by  a  purely  mathematical  argument, 
that  the  dependence  of  ckj  on  the  deformation  gradient 
history  must  be  of  the  following  form 
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t 

aij  =  xi,P  xj,Q  FPQ  [EABCt)]  » 

X=  -oo 

where  E^ft)  is  defined  by 


(2.8) 


£ABCt)  =  xk,A(T)xk,B(x)  '  6AB  »  (2*9) 

and  the  abbreviation  x.  p  is  used  for  the  deformation 
gradients  x^  p(t)  at  time  t.  The  nine  quantities  Eab(t) 
are  called  the  components  of  the  reduced  Cauchy  strain  ten¬ 
sor  at  time  t  and  we  note  that  EAB(t)  “  EBA(t)<  Accord¬ 
ingly,  only  six  of  the  quantities  Eab(t)  are  independent. 
Thus,  in  passing  from  (2.6)  to  (2.8),  we  have  replaced  arbit¬ 
rary  dependence  of  a..  on  the  nine  functions  x_  A(x)  by 
arbitrary  dependence  on  the  six  independent  components  of  the 
reduced  Cauchy  strain. 

So  far  we  have  made  no  assumption  regarding  the  symmetry 
of  the  material.  We  recall  that  in  classical  elasticity,  if 
the  mate  ial  possesses  some  symmetry,  restrictions  can  be 
placed  on  the  manner  in  which  the  stress  depends  on  the 
infinitesimal  strain  components.  In  the  case  when  the  material 
is  isotropic  these  restrictions  become  particularly  strong. 

In  principle,  we  can,  by  similar  considerations,  place 
restrictions  on  the  form  of  the  tensor-valued  functional 
FpQ  in  (2.8)  if  the  non-linear  viscoelastic  material  has  some 
symmetry.  It  emerges  that  the  functional  Fp^  must  satisfy 
the  condition 

t  t 

FPQ^\B^t^  =  aPMaQN  FMNrEAB^T)^  ’ 

f  —  -oo  T  =  *  00 


(2.10) 


9 


where 


EAB(t)  “  aAMaBNEMN(T)  * 


(2.11) 


for  all  belonging  to  the  group  of  transformations 

describing  the  symmetry  of  the  material.  In  the  case  when 
the  material  is  inotropic,  this  group  is  the  rotation  group, 
so  that,  for  an  isotropic  viscoelastic  material,  the  tensor¬ 
valued  functional  F^g  must  satisfy  the  condition  (2.10) 
for  all  rotations  a^  .  This  implicit  restriction  on  the 
form  of  FAB  can  be  made  explicit  for  any  specified  material 
symmetry  by  a  method  developed  by  Green  and  Rivlin  (1957)  and 
Wineman  and  Pipkin  (1964).  However,  we  shall  not  pursue  this 
here. 


If  the  material  considered  is  incompressible,  specifica¬ 
tion  of  the  deformation  gradient  history  does  not  determine 
the  stress  completely,  since  the  superposition  on  any  force 
system  of  a  hydrostatic  pressure  does  not  alter  the  deforma¬ 
tion.  Accordingly,  the  constitutive  assumption  (2.6)  must  be 
replaced  by 


:  VVaWJ 

X—  -oo 


(2.12) 


where  p  is  undetermined  if  the  deformation  gradient 
history  is  specified.  From  this  constitutive  assumption 
the  constitutive  equation 


t 

xi,Pxj,Q  FPQfEAB(x^ 

X  =  -  oo 


(2.13) 


follows,  replacing  (2.8,,  At  the  same  time,  the  fact  that, 


1JJ.W1L .  ,1 II I JJjppWflWg 


1^*4  W?PlWtns'-^ 


jin  ..iigg^n.Bjnn  . 
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in  an  incompressible  material,  the  volume  of  a  material 
element  must  remain  unchanged  during  deformation,  imposes  on 
the  reduced  Cauchy  strain  the  restriction 

det  lE^Cx)  +  5^1  =  1  .  (2.14) 

There  is  one  further  restriction  which  can  be  placed  on 
the  constitutive  equation  (2.8)  or  (2.13)  from  considerations 
of  a  very  general  character.  This  involves  the  concept  of  a 
hereditary  material  as  one  which  remains  unchanged  in  proper¬ 
ties  so  long  as  it  rests  undeformed.  For  such  a  material  a 
shift  in  time,  by  an  amount  tg  ,  say,  of  the  deformation 
history,  will  shift  the  dependence  of  stress  vs.  time  curve 
by  an  equal  amount  tg  parallel  to  the  time  axis. 

In  the  remainder  of  this  paper,  we  shall  restrict  our¬ 
selves  to  hereditary  materials. 
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3.  Non-linear  constitutive  equations 

In  §2  we  have  obtained  certain  restrictions  which  Bust 
be  inposed  on  the  manner  in  which  the  Cauchy  stress  depends 
on  the  deformation  gradient  hirtory.  These  are  expressed 
by  equation  (2.8)  in  which,  if  the  material  possesses  some 
symmetry,  the  tensor-valued  functional  Fpq  must  satisfy  the 
equation  (2.10). 

In  order  to  obtain  more  explicit  forms  for  the  constitu¬ 
tive  equation,  which  reflect  physically  non-pathological 
behavior  on  the  part  of  the  material  described,  various  rep¬ 
resentations  of  the  functional  Fpq  in  terms  of  series  of 
multiple  integrals  have  been  developed.  Green  and  Rivlin 
(1957,  1960)  and  Green,  Rivlin  and  Spencer  (1959)  considered 
materials  in  which  infinitesimal  changes  in  the  deformation 
history  result  in  infinitesimal  changes  in  the  stress  at 
time  t.  They  showed  that  a  sequence  of  approximations  to 
F pq  in  the  form  of  series  of  multiple  integrals  can  be  con¬ 
structed  which  tend,  in  the  limit,  to  Fpq.  In  order  to 
highlight  the  physical  content  of  this  development,  free  from 
the  mathematical  complications  which  arise  from  the  tensorial 
character  of  the  constitutive  equations ,  we  will  discuss  it 
in  the  context  of  time -dependent  simple  extension  of  a  rod. 

Accordingly,  let  o  be  the  stress  at  time  t  which  res¬ 
ults  from  a  time-dependent  fractional  extension  e(x).  Then, 
since  o  is  determined  if  e(x)  is  specified  for  --»<X£t  , 
we  may  say  that  a  is  a  functional  of  e(x),  thus: 
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t 

0  =  F  fe(T)]  . 

T=-<* 


(3.1) 


It  is  convenient  for  the  purposes  of  our  discussion  to  re¬ 
place  the  time  t  by  a  tine- like  variable  s  defined  by 


(t-x+l)p 


(3.2) 


where  p  is  positive.  Then,  for  a  specified  t,  the 

* 

extension  may  be  regarded  as  a  function  of  s  rather  than 
of  t  .  Noting  that  when  t  =  -  »  ,  s  =  0  and  when  x  =  t, 

s  =  1  ,  we  may  re-write  (3.1),  for  a  specified  t,  in  the 
form 


a 


1 

F 

s=0 


[e(s);t]  , 


(3.3) 


o  being  a  functional  of  e(s)  and  an  ordinary  function  of 
t.  For  hereditary  materials,  we  may  omit  the  dependence  of 
o  on  t  and  we  then  have 

1 

,a  ~  F  (e (s) ]  .  (3.4) 

s=0 

Now,  suppose  that  e(s)  can  be  represented  by  a  Fourier 
series  (see  §5)  thus 


e(s)  =  E  A  cos  n  Us  , 
n=0  n 


(3.5) 
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where 


n 


1 

e(s)cos  nils  ds  , 


0 


e(s)ds 


9 


0 


(3.6) 


The  condition,  under  which  the  representation  (3.5)  for 
e(s)  is  valid,  is  that  e(s)  have  bounded  variation  in 
the  range  s  =  0,  1  .  If  this  condition  is  satisfied,  then 
the  sum  of  the  series  (3.5)  gives  the  value  of  e(s)  at  all 
points  of  the  interval  s  =  0,  1  except  possibly  at  points 
of  discontinuity.  Accordingly,  if  we  assume  that  e(s)  is 
a  continuous  function  of  s,  then,  if  the  coefficients 
Aq,  Aj,.  . .  in  *iie  Fourier  series  (3.5)  are  specified,  the 
value  of  e(s)  is  determined  at  all  points  of  the  interval 
s  =  0,  1  and  o  may  be  regarded  as  a  function  of  Aq,A^,..., 
rather  than  as  a  function  of  e(s). 

We  now  suppose  that  in  (3.4),  a  is  a  continuous 

* 

functional  of  e(s)  ;  i.e.,  an  infinitesimal  change 


*  More  precisely,  if  e1(s)  and  e2(s)  are  two  extension 
histories,  and  and  a2  are  the  corresponding  values 

of  stress  at  time  t,  then  for  e  >  0  ,  there  exists  a 
value  of  6  >  0  ,  such  that 

ivo2i  <  c 

provided  that 

max  |e1  (s)  -  e2  (s) |  <6  . 
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in  e(s)  results  in  an  infinites imal  change  in  o  .  Then, 
a  nay  be  regarded  as  a  continuous  function  of  Aq,  A^,... 

To  some  degree  of  approximation,  we  cay  regard  o  as  a 
continuous  function  of  the  first  N  +  1  Fourier  coefficients 
Aq,  Aj,...,  A^  .  This  approximation  may  be  made  as  close  as 
we  please  by  taking  N  sufficiently  large.  We  denote  by 
such  an  approximation  to  o  .  Now,  from  Weiers trass 's 
theorem,  it  follows  that  we  can  approximate  as  closely 

as  we  please  by  a  polynomial  in  Aq,  A^,...,  A^  .  According¬ 
ly,  we  can  approximate  a  with  any  desired  accuracy  by  a 
polynomial  in  Aq,  A^,...,  Ajj  .  We  note  from  (3.6)  that  the 
product  of  r  ,  say,  A’s  is  an  t- tuple  integral  and  that  the 
kernel  in  this  integral  is  a  continuous  function  of  its  argu¬ 
ments.  For  example  , 


V2  ■ 


■IT 


COS  ns  1  cos 


2ns, 


e(s1)e(s2)ds1ds2  .  (3.7) 


Accordingly,  o  may  be  approximated  to  any  desired  accuracy 
by  the  sum  of  a  number  of  multiple  integrals,  thus 


*  fy(s1,s2,...,sy)e(s1)e(s2)...e(sy) 

0  0  _ 

ds^Sg. .  .dsy  =  o  (say) ,  (3.8) 


where  y  =  0,  1,  2,...  ,  and  the  kernels  in  the  multiple 
integrals  are  continuous  functions  of  s^.,  s2,...,  s  . 
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Hie  ten  in  (3.8)f  which  corresponds  to  u  =  0  ,  is  a 
constant.  This  nay  be  taken  to  be  zero,  if  we  assuae  that 
the  stress  is  zero  for  zero  strain  history  (i.e.,o  =0  when 
e(s)*0).  The  result  (3.8)  is  that  obtained  fro*  the  theory 
of  Green  and  Rivlin  (1957),  by  specializing  it  to  the  case 
of  simple  extension  of  a  rod. 


Using  (3.2),  we  can  rewrite  the  expressions  (3.6)  for 


A„  as 
n 


_a> 

--/• 


P  f  %  nil 

- ”+r  e(T)  cos  - o 

(t*T+l)P  1  (t-T+l)P 


dr  , 


(3.9) 


( t”  T+ 1) 


e(x)dT 


Correspondingly,  the  approximate  expression  for  a  in 
(3.8)  may  be  rewritten  as 


,T2,.  .  •  »TlJ) 


(3.10) 


e(x1)e(T2) . . .e(Ty)dT1dT2. . .dTy  , 


where 


^  p  (Ti  >t2  » •  •  • »  tjj) 


P^fy(Sl*S2>  *  *  *  ,SlP _ 

[(t-T1+l)  (t-T2+l)  .  .  .  (t-T^+1)]1 


(3.11) 


Since  fy(s1>s2 » • • • >sy)  is  a  bounded  function  of 
slts2  ...,s^  ,  and  hence  of  t1  ,t2>...,t^  ,  it  follows 
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that  fjCr  ,t  becomes  vanishingly  snail  as  any  of 

the  variables  i^,  t2,...,tv  tends  to  -®.  This  inplies 
that  the  naterial  has  fading  nenory  in  the  usual  sense  that, 
as  t  -  t  +  ®  ,  the  effect  on  a  of  the  extension,  in  a 
finite  tine  interval  about  tine  t  ,  becomes  vanishingly  snail. 

Fron  a  purely  nathenatical  point  of  view,  if  e(s) 
possesses  one  or  nore  discontinuities  in  the  interval  s  =  0,1, 
then  in  order  to  determine  e(s)  at  all  points  of  this  inter¬ 
val,  we  nust  specify,  not  only  the  values  of  tho  Fourier  co¬ 
efficients  given  by  (3.6),  but  also  the  values  of  e(s)  at 
the  various  values  of  s  at  which  the  discontinuities  occur. 

Let  us  suppose  that  discontinuities  occur  at  s  *  S^(A=1,2,. . .  ,v) 
and  let  e^=  e(S^)  •  Then,  if  o  is  a  continuous  functional 
of  e(s)  ,  we  may  regard  it  as  a  continuous  function  of 
Aq,  A^,...  and  of  e^ (A=l,2 , . . .  ,v)  ,  the  nature  of  this 
function  depending  on  the  values  of  s  at  which  the  dis¬ 

continuities  occur.  Approximate  expressions  o  for  a  may 
then  bv  obtained  in  the  form  (3.8)  ,  where  f  is  also  a 
function  of  e ^  and  ,  as  well  as  of  s1,s2,..,,s  . 

The  term  corresponding  to  y  =  0  is  now  a  function  of 
e^  and  , 

In  considering  the  applicability  of  this  conclusion  to 
real  materials,  it  is  well  to  bear  in  mind  that  discontinuous 
changes  in  the  extension  cannot  in  fact  be  produced.  The 
importance  of  considering  a  discontinuous  change  in  e(s)  at 
time  5  =  5.^,  say,  resides  only  in  the  fact  that  it  provides 
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an  idealized  model  for  a  very  rapid,  but  continuous,  change 
in  the  value  of  e(s)  in  the  interval  S1~  e  ,  S.+  e  , 
where  e  is  snail.  Consequently,  we  nay,  in  general,  onit 
consideration  of  discontinuities  in  e(s).  An  exception 
arises  when  this  discontinuity  occurs  at  the  instant  at  which 
the  stress  is  neasured,  i.e.,  at  s  =  1  .  In  order  to  accomno- 
date  perfectly  elastic  materials  and  materials  exhibiting 
instantaneous  elasticity  within  the  franework  of  the  theory, 
we  nust  then  include,  in  the  expression  for  o  ,  explicit 
dependence  on  the  instantaneous  value  e  of  e(s)  at  s  *  1. 
Approximate  expressions  o  for  a  nay  then  be  obtained  in 
the  form  (3.8),  where  fy  is  a  function  of  e  ,  as  well  as  of 
sx,  s2,...fs^  .  The  term  corresponding  to  p  =  0  is  now  a 
function  of  e  .  This  result  is  essentially  that  obtained  from 
the  theory  of  Green,  Rivlin  and  Spencer  (1959),  by  specializ¬ 
ing  it  to  the  case  of  simple  extension  of  a  rod. 

It  is  perhaps  worth-while  to  underline  the  fact  that, 
whether  or  not  the  f's  depend  on  e  ,  the  approximation  to 
a  represented  by  (3.8)  is  of  the  following  type.  If  e  is 
any  specified  positive  quantity,  we  can  construct  an  approxi¬ 
mation  o  ,  say,  to  o  of  the  form  (3.8),  such  that 

\a  -  a j  <  e  (3.12) 

for  all  extension  histories  e(s)  which  have  bounded 


variation. 


The  kernels  occuring  in  the  expression  for  a  cannot 
be  uniquely  determined.  We  can,  however,  construct  a  sequence 
of  approximations  of  the  form  (3. 8) ,  corresponding  to  smaller 
and  smaller  c  ,  which  tends  in  the  limit  tc  o  .  Corres¬ 
ponding  kernels  in  this  sequence  do  not  necessarily  tend  to  a 
limit  so  that  we  cannot  say  that  a  may,  without  error,  be 
represented  by  an  expression  of  a  form  similar  to  that  in 
(3.8).  The  absolute  error  involved  in  approximating  o  by 
a  may,  however,  be  made  as  smai!  as  we  please.  In  the  neigh¬ 
borhood  of  a  =  0  ,  the  percentage  error  may  then  be  large. 

It  is,  however,  perfectly  safe  to  use  a  representation 
for  o  of  the  form  (3. 8), provided  that  we  are  concerned  only 
with  changes  of  0  of  magnitude  much  greater  than  e  .  Since 
e  may  be  made  as  small  as  we  please,  it  might,  at  first  sight, 
appear  that  this  does  not  impose  a  meaningful  restriction  on 
the  use  of  the  representation  (3.8).  Let  us  therefore  consid¬ 
er  an  example  in  which  this  might,  in  fact,  impose  a  meaning¬ 
ful  restriction.  Suppose  we  consider  extension  histories 
ae(s)  ,  with  e  zero,  where  a  may  be  made  as  small  as  we 
please.  For  small  enough  a  ,  and  taking  fg=  0  ,  we  may 
approximate  a  by  the  fi^st  term  in  the  series  (3.8)  thus 

1 

f^s^eCs^dSj^  .  (3.13) 

This  approximation  becomes  increasingly  good  as  as  a 
decreases.  However,  regarded  as  an  approximation  to  o  , 
the  expression  on  the  right-hand  side  of  (3.13)  has  th? 
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limitation  that  a  only  approximates  a  with  absolute  error 
less  than  e  .  It  nay  be  that,  in  order  for  (3.13)  to  pro¬ 
vide  a  good  approximation  to  a  ,  the  value  of  a  must  be  so 
decreased  that  the  stress  is  comparable  with  c  .  The  fact 
that  e  may  be  made  as  small  as  we  please  does  not  save  us 
from  this  limitation,  since  choice  of  a  lower  value  of  c 
requires  a  new  approximate  representation  of  the  form  (3.8) 
and  with  this  we  nay  be  driven  to  lower  values  of  a  and 
hence  of  the  stress. 

If,  however,  a  may  be  expressed  exactly  in  the  fora 
(3.8),  then,  for  sufficiently  small  values  of  a  ,  we  may 
approximate  a  by  the  expression  (3.13)  and,  for  somewhat 
larger  values,  by 


-  « f1  f1(s1)e(s1)ds1  *  jf1  f1f2(s1,s2) 
j  f\  J  i\  J{\ 


0  •'0 

e(s1)e(s2)ds1ds2  ,  (3.14) 


and  so  on.  In  this  way,  we  obtain  a  hierarchy  of  expressions 
for  a  ,  valid  over  larger  and  larger  ranges  of  the  extension 
prior  to  time  t  .  This  is,  of  course,  also  true  if  a  can 
be  approximated  by  an  expression  of  the  form  (3.8),  with  an 
error  which  is  of  degree  greater  than  p  ,  say,  in  the  exten¬ 
sion.  Of  course,  whether  or  not  one  of  the  approximations 
has  significant  value  in  the  study  of  a  particular  material, 
to  which  it  applies  in  principle,  will  depend  on  whether  the 

range  of  extensions  for  which  j t  is  valid  is  a  range  over 
which  experiments  can  be  carried  out  with  reasonable  accuracy. 
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Coleman  sad  Noll  (I960)  hove  obtained  n  successive 
approximations  essentially  of  the  forms  (3.13)  and  (3.14) 
by  assuming  that  the  fractional  F  in  (3.3)  is  Frdchet 
differentiable  n  times  about  the  zero  history  (i.e. ,  the 
history  e(s)  *  0  ,  0<s£l) .  Their  argument  is  basically 

circular  in  that  the  definition  of  nth  order  Frddtet 
differentiability  is  that  the  functional  F  shall  be  ex¬ 
pressible  as  the  sum  of  homogeneous  functionals  of  degrees 

1,  2 . n  in  e(s)  ,  with  a  residue  of  degree  greater 

than  n  in  e(s)  . 

The  discussion  has*  so  far,  been  based  on  the  assumption 
that  the  stress  at  tine  t  is  a  continuous  functional  of  the 
extension  history  e(x)  for  -  •»  <  x£t  .  It  is  evident  that 
this  is  not  necessarily  an  appropriate  assunption  for  many  of 
the  materials  with  which  we  are  concerned.  For  example,  let 
us  consider  two  deformations  in  both  of  which  the  extension 
is  maintained  at  zero  up  to  and  including  time  t  .  In  one 
of  these,  we  maintain  the  extension  zero  after  time  t  ,  while 
in  the  other  the  rate  of  extension  is  changed  dis continuously 
at  time  t  from  zero  to  a  finite  value.  In  many  materials 
(e.g. ,  in  Newtonian  fluids)  the  stress  at  time  t  will  be 
quite  different  in  the  two  cases. 

In  order  to  accommodate  such  materials  in  the  framework 
of  our  mathematics,  we  may  take  our  constitutive  assumption 
in  the  following  form  :  the  stress  is  a  functional  of  e(t) 
up  to  and  including  time  t  and  an  ordinary  function  of 
e  (the  rate  of  extension  at  time  t  )  ,  thus  : 

a  -  V  [e(x)  ;  e]  .  (3.15) 


Generalizing  this  concept  to  accommodate  materials  for  which 
discontinuous  changes  at  time  t  in  the  higher  time  deriva¬ 
tives  e  ,  e,...  of  the  extension  affect  the  stress  at  time 
t  ,  we  may  make  our  initial  constitutive  assumption  in  the 
form  :  the  stress  at  time  t  is  a  functional  of  e(x)  in 
the  interval  -  «  <  t  <  t  and  an  ordinary  function  of 
e  ,  e , . . . ,  thus 

t 

a  =  F  [e(x)  ;  e,  e,...  ]  .  (3.16) 

~oo 

This  is  esentially  the  constitutive  assumption  used  by 
Green  and  Rivlin  (1960),  specialized  to  the  case  of  simple 
extension  of  a  rod. 

Alternatively,  we  could  accommodate  such  constitutive 
assumptions  in  the  single  assumption  that  the  stress  a  at 
time  t  is  given  by 

t  t+e 

o  =  Lz  F  [e(x)J  ,  (3.17) 

£->0  T=-oo 

for  small  positive  e  ,  i.e.,  stress  at  time  t  is  the 
limit  as  e  0  of  a  functional  of  e(i)  over  the  range 
-oo<x<t+e.  To  see  this  we  note  that 


e  =  -  L 
e-+0 


e  = 


L 

€-*■0 


/t+e 

-00 

t  rt+£ 

/  6" (t- t)e  (r 

°J 


(T-t)e(t)dr  , 

) di  ,  etc. , 


(3.18) 


where  6(  )  denotes  the  Dirac  delta  function  and  S'(  ) 
denotes  its  derivative.  It  should  be  noted  that  even  if, 
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in  (3.15)  and  (3.16),  a  is  a  continuous  functional  of 
o(t)  and  a  continuous  function  of  the  instantaneous  values 
of  the  time  derivatives  of  e  at  time  t  ,  the  functional 
dependence  on  e(t)  expressed  by  (3.17)  is  not,  in  general, 
continuous . 

Of  course,  if  we  limit  ourselves  to  deformations  in 
which  the  time  derivatives  of  e(t)  are  continuous  at  t  , 
then  the  constitutive  assumptions  (3.15)  and  (3.16)  can  be 
replaced  by  (3.1). 

Approximations  to  a  in  the  form  of  series  of  multiple 
integrals  may  be  made  in  the  cases  when  o  is  given  by 
(3.15)  or  (3.16).  The  approximation  takes  the  form  (3.10), 
with  the  kernels  dependent  on  e  in  the  case  (3.15)  and  on 
e  ,  e,...  in  the  case  (3.16). 

It  should  be  appreciated,  in  interpreting  the  above 
remarks,  that  if  the  '’slope”  of  e(r)  vs.  t  changes  dis- 
continuously  at  time  t  ,  then  strictly  its  time  derivative 
at  time  t  does  not  exist.  B/  e  we  mean  the  right-hand 
derivative  of  e(x)  at  time  t  .  This  is  in  accord  with 
the  usual  convention  in  mechanics  when  we  apply  the  Navier- 
Stokes  equation,  say,  to  situations  in  which  the  rate  of 
deformation  changes  discontinuously  at  time  t  . 


I 
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4.  Generalization  to  arbitrary  deformations 

In  the  previous  section,  we  have  attenpted  to  show 
how  mathematical  expression  nay  be  given  to  a  variety  of 
possible  types  of  mechanical  behavior  which  nay  be  exhibited 
by  viscoelastic  materials.  We  have  done  this  in  the  rela¬ 
tively  simple  context  of  simple  extension. 

FOr  each  type  of  behavior  we  may  obtain  analogous 
equations  which  are  valid  for  deformations  of  a  general 
character.  Thus,  if  we  take  our  initial  constitutive  assump¬ 
tion  in  the  form  (2.6),  we  have  seen  that  the  Cauchy  stress 
must  necessarily  be  expressible  in  the  form  (2.8).  Replacing 
the  time  r  by  the  time -like  variable  s  ,  defined  by 
(3.2),  we  see  that  the  stress  ,  referred  to  a  rectangu¬ 

lar  Cartesian  system  x  ,  must  be  expressible  in  the  form 

1 

aij  =  Xi,PXj,Q  FPQ  [EPQ(s)1  *  (4,1) 

s=0 


where  Epq(s)  is  the  reduced  Cauchy  stress  at  time  s  . 

This  can  be  re-written  more  succinctly  in  matrix  notation, 
thus  : 

1  y 

a  =  F  F  [ECs) JF1  ,  (4.2) 

s=0 


where 


a 


n°ij 


ECs)  =  I  lE^fs)  |  | 


(4.3) 
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T 

and  F  denotes  the  transpose  of  F  . 

We  now  suppose  that  E(s)  can  be  represented  by  a 
Fourier  series  thus 


E(s)  =  An  cos  nils  ,  (4.4) 


where  An  is  given  by 


E(s)ds  ; 


(4.5) 


i.e.,  each  of  the  components  of  the  tensor  E(s)  can  be 
expressed  as  a  Fourier  series.  This  is  possible  under  the 
same  conditions  as  apply  in  our  discussion  of  the  case  of 
simple  extension.  Each  of  the  coraponc its  of  E(s)  must 
have  bounded  variation  in  the  interval  [0,1].  Paralleling 
the  discussion  in  §2,  we  assume  that  each  of  the  compon¬ 
ents  of  E(s)  is  continuous  everywhere,  except  possibly  at 
s  =  1  .  Then,  in  order  to  represent  E(s)  over  the 
complete  range  s  =  [0,1]  ,  we  must  specify  not  only  the 
coefficients  An  in  (4.4),  but  also  the  value  E  of  E(s) 
at  s  =  1. 


The  tensor  F  in  (4.2)  is,  of  course,  a  tensor-valued 
functional  and  we  suppose  that  it  is  such  that  infinitesimal 
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changes  in  the  components  of  E(s)  lead  only  to  infinitesi¬ 
mal  changes  in  the  components  of  F  .  The  tensor-valued 
functional  F  is  then  said  to  be  a  continuous  functional  of 
the  tensor  E(s)  .  Analogously  with  the  case  of  simple  exten¬ 
sion,  we  may  approximate  each  of  the  components  of  F  to  any 
desired  accuracy  by  a  polynomial  in  the  components  of  the 
tensors  E  and  .  Thus,  each  of  the  components  Fp^  of 
F  may  be  expressed,  with  any  desired  accuracy,  in  the  form 


l  “pQAiV 


.A  B 

y  y 


<V  ("2> 

Vl  y»2-- 


A°V 

V* 


9 


“  TpQ  (say) 


(4.6) 


where,  for  specified  n,  are  the  components  of 

A  ;  also  n,  ,  n„,...  take  integral  values  and  the  a's 

„n  12 

are  functions  (or,  if  we  like,  polynomial  functions)  of  the 
components  of  E  . 

Using  (4.5),  we  may  rewrite  (4.6)  in  the  form 


F  ^  F 
r:q  ~  rPQ 


=  £  J  *••  /  fPQA,B1...AB  (si»--*sy) 

V  Jn  Jn  11  ^ 


EA1B1(s1)  •  •  ’  ,dsy  ’ 


(4.7) 


where  the  kernels  fDnA  r  A  R  are  continuous  functions 

PQAlBl-'’Vy 

of  s1,...,sy  and  functions  (or  polynomial  functions)  of 
the  components  E^g  of  E  . 
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Using  (3.2),  we  can,  of  course,  re-write  (4.7)  as 


FPQ  '  fPQ  "  £ 


J  -J  TPQA1B1.  -A^B^l . V 


EAiB1^Tl^***EAiiBi|^Ty^dTl*"dTp  *  (4*8) 


where 


“%QA1B1...A11B1| 


(4.9) 


If  the  material  considered  has  some  symmetry,  the 
functional  F  must  satisfy  the  condition  (2.10).  In 
approximating  F  by  T  ,  to  any  desired  accuracy,  this 
can  be  done  by  an  expression  of  the  form  (4.7)  which  satis¬ 
fies  the  restriction  (2.10)  imposed  by  material  symmetry. 

We  can  develop  constitutive  equations  for  arbitrary 
deformations  corresponding  to  the  simple  extensional  case 
presented  in  (3.16).  We  start  with  the  constitutive  assump¬ 
tion  that  tiie  Cauchy  stress  components  at  time  t  are 
functionals  of  the  histories  of  the  deformation  gradients 
x  a(t)  f°r  -  00  <  t  <  t  and  functions  of  the  velocity 

p  ,/v  — 

gradients  Xp  ^  ,  acceleration  gradients  3^  ^  ,  and  so  on,- 
at  time  t  .  Thus , 


a. .  =  F. . [x  ,(x) 
ij  i  i  1  p ,  A 


P,A  ’  p , A 


(4.10) 
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Paralleling  the  passage  from  (2.6)  to  (2.&),  we  reach  the 
conclusion,  in  this  case,  that  the  stress  nust  be  express¬ 
ible  in  the  fora 


fij  '  xi,Pxj,QFPQfE 


AB 


(t) 


^AB 

~at" 


d2E 


AB 


dt* 


(4.11) 


Using  the  notation  (4.3),  we  can,  of  course,  rewrite 
(4.10)  as 

dE  d2E 

o  =  FF[E(s)  J  ST  ,  — Z  •  (4.12) 

at  dt2 

If  the  dependence  of  Fpq  on  E^^)  is  continuous, 
then  we  can  approximate  Fpq  with  any  desired  accuracy  by 
a  series  of  multiple  integrals  in  the  form  (4.8),  the  ker¬ 
nels  now  being  functions  of  the  components  of  the  reduced 
Cauchy  strain  and  its  time  derivatives  at  time  t  .  If  the 
dependence  of  Fpq  on  E  ,  dEAg/dt  ,  d2EAfi/dt<1 , . . .  is 
continuous,  then  the  dependence  of  the  kernels  on  these  may 
be  taken  as  polynomial,  with  any  desired  accuracy. 

Just  as  in  the  previous  cases  discussed,  if  the  mater¬ 
ial  considered  has  some  symmetry,  the  tensor- valued 
functional  in  (4.10)  must  satisfy  the  restrictions  implied 
by  (2.10). 


In  order  to  arrive  at  the  Fourier  representation 
(3.5)  for  e(s)  ,  we  proceed  in  the  following  Banner. 
Me  first  fora  a  function  e(s)  given  by 


e(s)  =  e(s)  s  =  [0,1]  , 

e(s)  *  e(-s)  s  =  (-1,0)  . 


(5.1) 


e(s)  is  thus  an  even  function  of  s  ,  defined  in  the 

e 

range  [-1,1]  by  (5.1). 

We  now  define  e(s)  ,  not  just  for  s  =  [-1,1]  , 
but  for  all  values  of  s  ,  as  a  periodic  function  with 
periodicity  2,  which  is  given  by  (5.1)  in  the  range 
[-1.1]  .  We  express  this  as  a  cosine  Fourier  series  thus. 


e(s)  =  l  A  cos  nIIs  » 
n=0  n 


(5.2) 


where,  using  the  first  of  equations  (5.1)  , 


f1  -  r1 

A n  =  1  e(s)  cos  nils  ds  =  2  I 

•'-l  Jo 

Aq  =  jJ  e(s)ds  =  J  e(s)ds 


e(s)  cos  nils  ds  =  2  J  e(s)  cos  nils  ds  (n>^l)  , 

Jo 

(5.3) 


In  view  of  the  first  of  equations  (5.1),  we  obtain  (3.5)  . 


This  paper  was  written  with  the  support  of  the  Office 
of  Naval  Research  under  Contract  No.  N00014-67-0370-0001 
with  Lehigh  University.  My  thanks  are  due  to  Or.  D.G.B. 
Edelen  for  interesting  discussions  in  connection  with  this 
work. 
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